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We survey the critical behaviour of strongly interacting matter, as determined in 
finite temperature lattice QCD. 

Keywords: 05.70.Fh; 12.38.Mh. 

1. Phases of Strongly Interacting Matter 

What happens to strongly interacting matter in the limit of high temperature 
and/or density? This question has fascinated physicists ever since short range strong 
interactions and the resulting multiple hadron production were discovered. Let us 
look at some of the features that have emerged. 

• Hadrons have an intrinsic size, with a radius — I fm, and hence a hadron 
needs a space of volume Vh — (47r/3)r?in order to exist. This suggests a 
limiting density n c of hadronic matteff, with n c — 1/Vh — 1.5 no, where 
no ~ 0.17 fm~ 3 denotes the density of normal nuclear matter. 

• Hadronic interactions provide abundant resonance production, and the re- 
sulting number p(m) of hadron species increases exponentially as function 
of the resonance mass m, p(m) ~ exp(6m). Such a form for p(m) appeared 
first in the statistical bootstrap model, based on self-similar resonance for- 
mation or decajEl It was then also obtained in the more dynamical dual 
resonance approach^. In hadron thermodynamics, the exponential increase 
of the resonance degeneracy results in an upper limit for the temperature 
of hadronic mattell, T c = 1/6 ~ 150 - 200 MeV. 

• What happens beyond T c l In QCD, hadrons are dimensionful color-neutral 
bound states of more basic pointlike colored quarks and gluons. Hadronic 
matter, consisting of colorless constituents of hadronic dimensions, can 
therefore turn at high temperatures and/or densities into a quark-gluon 
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plasma of pointlike colored quarks and gluons as constituent^. This dc- 
confincment transition leads to a colour-conducting state and thus is the 
QCD counterpart of the insulator-conductor transition in atomic matted 

• A shift in the effective constituent mass is a second transition phenomenon 
expected from the behavior of atomic matter. At T = 0, in vacuum, quarks 
dress themselves with gluons to form the constituent quarks that make 
up hadrons. As a result, the bare quark mass m q ~ is replaced by a 
constituent quark mass M q ~ 300 MeV. In a hot medium, this dressing 
melts and M q — > 0. Since the QCD Lagrangian for m q — is chirally 
symmetric, M q ^ implies spontaneous chiral symmetry breaking. The 
melting M q — > thus corresponds to chiral symmetry restoration. We shall 
see later on that in QCD, as in atomic physics, the shift of the constituent 
mass coincides with the onset of conductivity. 

• A third type of transition would set in if the attractive interaction between 
quarks leads in the deconfined phase to the formation of colored bosonic 
diquark pairs, the Cooper pairs of QCD. These diquarks can then condense 
at low temperature to form a color superconductor. Heating will dissociate 
the diquark pairs and turn the color superconductor into a normal color 
conductor. 

Using the baryochemical potential \i as a measure for the baryon density of the 
system, we thus expect the phase diagram of QCD to have the schematic form shown 
in Fig. ^ Given QCD as the fundamental theory of strong interactions, we can use 
the QCD Lagrangian as dynamics input to derive the resulting thermodynamics of 
strongly interacting matter. For vanishing baryochemical potential, fj, = 0, this can 
be evaluated with the help of the lattice regularisation, leading to finite temperature 
lattice QCD. 



T 




hadronic matter 



quark-gluon plasma 



diquark matter 



Fig. 1. The phase diagram of QCD 
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2. Finite Temperature Lattice QCD 

QCD is defined by the Lagrangian 

£ = - \F% F a U - E tKrfd* + mj- fffA,,)^ , (1) 
/ 

with 

F% = (d,Al - duAl - gK c A\K) ■ (2) 

Here A a ^ denotes the gluon field of colour a (a=l,2,...,8) and ip^ the quark field of 
colour a (a=l,2,3) and flavour /; the input ('bare') quark masses are given by m/. 
With the dynamics thus determined, the corresponding thermodynamics is obtained 
from the partition function, which is most suitably expressed as a functional path 
integral, 

Z(T,V) = J dA dip drp exp (- J d 3 x ' dr C{A,ip,$)\ , (3) 

since this form involves directly the Lagrangian density defining the theory. The 
spatial integration in the exponent of Eq. Q is performed over the entire spatial 
volume V of the system; in the thermodynamic limit it becomes infinite. The time 
component xq is "rotated" to become purely imaginary, r = ixo, thus turning the 
Minkowski manifold, on which the fields A and tp are originally defined, into a 
Euclidean space. The integration over r in Eq. © runs over a finite slice whose 
thickness is determined by the temperature of the system. 

From Z(T, V), all thermodynamical observables can be calculated in the usual 
fashion. Thus 



gives the energy density, and 



dV 



T 



the pressure. For the study of critical behaviour, long range correlations and multi- 
particle interactions are of crucial importance; hence perturbation theory cannot 
be used. The necessary non-perturbative regularisation scheme is provided by the 
lattice formulation of QCEpl it leads to a form which can be evaluated numerically 
by computer simulatiorpl. 

The first variable to consider is the deconfincmcnt measure given by the Polyakov 
loop™ 

L(T) - lim exp{-Wr)/T} (6) 

r — >oo 

where V(r) is the potential between a static quark- antiquark pair separated by a 
distance r. In the limit of large input quark mass, V(oo) = oo in the confined phase, 
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so that then L = 0. Colour screening, on the other hand, makes V(r) finite at large 
r, so that in the deconfined phase, L does not vanish. It thus becomes an 'order 
parameter' like the magnetisation in the Ising model. In the large quark mass limit, 
QCD reduces to pure SU(3) gauge theory, which is invariant under a global Z 3 
symmetry. The Polyakov loop provides a measure of the state of the system under 
this symmetry: it vanishes for Z3 symmetric states and becomes finite when Z3 is 
spontaneously broken. Hence the critical behaviour of SU(3) gauge theory is in the 
same universality class as that of Z3 spin theory (the 3-state Potts model): both 
are due to the spontaneous symmetry breaking of a global Z3 symmetry^! 

For finite quark mass m q , V(r) remains finite for r — > 00, since the 'string' 
between the two colour charges 'breaks' when the corresponding potential energy 
becomes equal to the mass Mh of the lowest hadron; beyond this point, it becomes 
energetically more favourable to produce an additional hadron. Hence now L no 
longer vanishes in the confined phase, but only becomes exponentially small there, 



L(T) ~ exp{-M h /T}; 



(7) 



here Mh is of the order of the p-mass, so that L ~ 10~ 2 , rather than zero. Decon- 
finement is thus indeed much like the insulator-conductor transition, for which the 
order parameter, the conductivity cr(T), also does not really vanish for T > 0, but 
with <t(T) ~ exp{— AE/T} is only exponentially small, since thermal ionisation 
(with ionisation energy AE) produces a small number of unbound electrons even 
in the insulator phase. 

Fig. |5]a shows recent lattice result*^ for L(T) and the corresponding suscep- 
tibility Xl(T) ~ (L 2 ) — (L) 2 . The calculations were performed for the case of two 
flavours of light quarks, using a current quark mass about four times larger than 
that needed for the physical pion mass. We note that L(T) undergoes the expected 
sudden increase from a small confinement to a much larger deconfinement value. 
The sharp peak of xl{T) defines quite well the transition temperature T c , which 
we shall shortly specify in physical units. 
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Fig. 2. Polyakov loop and chiral condensate in two-flavour 
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The next quantity to consider is the effective quark mass; it is measured by the 
expectation value of the corresponding term in the Lagrangian, (ipip) (T) . In the limit 
of vanishing current quark mass, the Lagrangian becomes chirally symmetric and 
the corresponding order parameter. In the confined phase, with effective 
constituent quark masses M q ~ 0.3 GeV, this chiral symmetry is spontaneously 
broken, while in the deconfined phase, at high enough temperature, we expect its 
restoration. In the real world, with finite pion and hence finite current quark mass, 
this symmetry is also only approximate, since (ipip)(T) now never vanishes at finite 
T. 

The behaviour of ('tpifj) (T) and the corresponding susceptibility Xm ~ 
d^ip) /dm q are shown in Fig. 0b 13 , calculated for the same case as above in Fig. 
0a. We note here the expected sudden drop of the effective quark mass and the 
associated sharp peak in the susceptibility. The temperature at which this occurs 
coincides with that obtained through the deconfinement measure. We therefore con- 
clude that at vanishing baryon number density, quark deconfinement and the shift 
from constituent to current quark mass coincide. We thus obtain for ^ = a 
rather well defined phase structure, consisting of a confined phase for T < T c , with 
L(T) ~ and (iJip)(T) ^ 0, and a deconfined phase for T > T c with LIT) ^ and 
{ipil}){T) ~ 0. The underlying symmetries associated to the critical behaviour at 
T = T Cl the Z% symmetry of deconfinement and the chiral symmetry of the quark 
mass shift, become exact in the limits m q — > oo and m q —> 0, respectively. In the 
real world, both symmetries are only approximate; nevertheless, we see from Fig. [5] 
that both associated measures retain an almost critical behaviour. 

Next we come to the behaviour of energy density e and pressure P at deconfine- 
ment. In Fig. |2| it is seen that e/T 4 changes quite abruptly at the above determined 
critical temperature T = T c , increasing from a low hadronic valu e to one slightly 
below that expected for an ideal gas of massless quarks and gluon^3. 




Besides the sudden increase at deconfinement, there are two further points to 
note. In the region T c < T < 2 T c , there still remain strong interaction effects. As seen 
in Fig.^J the 'interaction measure' A = (e — 3P)/T 4 remains sizeable and does not 
vanish, as it would for an ideal gas of massless constituents^!. In the simple model 
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of the previous section, such an effect arose due to the bag pressure, and in actual 
QCD, one can also interpret it in such a fashion^. It has also been considered in 
terms of a gradual onset of deconfinement starting from high mome ntaj ^j and most 
recently as a possible consequence of coloured "resonance" states^. The second 
point to note is that the thermodynamic observables remain about 10% below their 
Stefan-Boltzmann values (marked "SB" in Fig. EJ) even at very high temperatures, 
where the interaction measure becomes very small. Such deviations from ideal gas 
behaviour can be expressed to a large extent in terms of effective 'thermal' masses 
m t h of quarks and gluons, with m t h ~ g(T) ylEEDEI] Maintaining the next-to- 
leading order term in mass in the Stefan-Boltzmann form gives for the pressure 



P = cT 4 

and for the energy density 
e = 3cT 4 

= 3cT 4 



3 V T ) 
l-ag 2 (T) + 



= cT*[l-ag 2 (T)} 



(8) 



2a / m tt 
3 V T 

2a mth ( dg 
dT 



( dm 



lh 



V dT 



(9) 



where c and a are colour- and flavour-dependent positive constants. The deviations 
of P and e from the massless Stefan-Boltzmann form thus vanish only logarithmi- 
cally, with g(T) ~ 1/ln T, while the interaction measure 

,dg 



A = 2acg{T) T 



dT 



(10) 



vanishes up to terms of order Tdg/dT. 

Finally we turn to the value of the transition temperature. Since QCD (in the 
limit of massless quarks) does not contain any dimensional parameters, T c can 
only be obtained in physical units by expressing it in terms of some other known 
observable which can also be calculated on the lattice, such as the p-mass, the 
proton mass, or the string tension. In the continuum limit, all different ways should 
lead to the same result. Within the present accuracy, they define the uncertainty 
so far still inherent in the lattice evaluation of QCD. Using the p-mass to fix the 
scale leads to T c ~ 0.15 GeV, while the string tension still allows values as large as 
T c ~ 0.20 GeV. This means that energy densities of some 1-2 GeV/fm 3 are needed 
in order to produce a medium of deconfined quarks and gluons. 

In summary, finite temperature lattice QCD at vanishing baryon density shows 

• that there is a transition leading to colour deconfinement coincident with 
chiral symmetry restoration at T c ~ 0.15 - 0.20 GeV; 

• that this transition is accompanied by a sudden increase in the energy 
density (the "latent heat of deconfinement") from a small hadronic value 
to a much larger value, about 10% below that of an ideal quark-gluon 
plasma. 
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In the following section, we want to address in more detail the nature of the critical 
behaviour encountered at the transition. 

3. The Nature of the Transition 

We begin with the behaviour for vanishing baryon density (fi — 0) and come to 
/i =/= at the end. Consider the case of three quark species, u, d, s. 

• In the limit m q — > oo for all quark species, we recover pure SU(3) gauge 
theory, with a deconfincment phase transition provided by spontaneous Z3 
breaking. It is first order, as is the case for the corresponding spin system, 
the 3-state Potts model. 

• For m q — > for all quark masses, the Lagrangian becomes chirally symmet- 
ric, so that we have a phase transition corresponding to chiral symmetry 
restoration. In the case of three massless quarks, the transition is also of 
first order. 




crossover region 








CO 



Fig. 5. The nature of thermal critical behaviour in QCD 

• For < m q < 00, there is neither spontaneous Z 3 breaking nor chiral sym- 
metry restoration. Hence in general, there is no singular behaviour, apart 
from the transient disappearence of the first order discontinuities on a line 
of second order transitions. Beyond this, there is no genuine phase transi- 
tion, but only a "rapid cross-over" from confinement to deconfincment. The 
overall behaviour is summarized in Fig. [S] 

• As already implicitely noted above, both "order parameters" L(T) and 
x(T) nevertheless show a sharp temperature variation for all values of m qi 
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so that it is in fact possible to define quite well a common cross-over point 
T . 

• The nature of the transition thus depends quite sensitively on the number of 
flavours Nf and the quark mass values: it can be a genuine phase transition 
(first order or continuous), or just a rapid cross-over. The case realized 
in nature, the "physical point" , corresponds to small u, d masses and a 
larger s-quark mass. It is fairly certain today that this point falls into the 
cross-over region. 

• Finally we consider briefly the case of finite baryon density, /i / ^2. 
Here the conventional computer algorithms of lattice QCD break down, and 
hence new calculation methods have to be developed. First such attempts 
(re weighting^l, analytic continuation^! power serieJ^l) suggest the phase 
diagram shown in Fig. El The position of the critical end-point of the line 
of first order transition depends on the position of the physical point in the 
mu,d ~ m s plane, as illustrated in Fig. [5] 



T 




Fig. 6. Deconfincmcnt in the T,fi plane 



The critical behaviour for strongly interacting matter at low or vanishing baryon 
density, describing the onset of confinement in the early universe and in high energy 
nuclear collisions, thus occurs in the rather enigmatic form of a "rapid cross-over" . 
There is no thermal singularity and hence, in a strict sense, there are neither dis- 
tinct states of matter nor phase transitions between them. So what does the often 
mentioned experimental search for a "new state of matter" really mean? How can a 
new state appear without a phase transition? Is there a more general way to define 
and distinguish different states of bulk media? After all, in statistical QCD one does 
find that thermodynamic observables - energy and entropy densities, pressure, as 
well as the "order parameters" L(T) and x(T) ~ continue to change rapidly and 
thus define a rather clear transition line in the entire cross-over region. Why is this 
so, what is the mechanism which causes such a transition? 
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In closing this section, we consider a speculative answer to this rather fundamen- 
tal question. It is known that in spin systems, when the partition function becomes 
analytic and hence there is no more thermal critical behaviour, the onset of large 
cluster formation, i.e., percolation, persists as a geometric form of critical behaviour 
persitJ^HI. This can be applied also to the question of deconfinemenl 29 . Consider 
hadrons of intrinsic size Vh = (47r/3)r^, with ~ 0.8 fm. In three-dimensional 
space, the formation of a connected large-scale cluster first occurs at the density 

n c = ^ ~ 0.16 fm" 3 . (11) 
Vh 

This point specifies the onset of hadronic matter, in contrast to a gas of hadrons, 
and it indeed correctly reproduces the density of normal nuclear matter. However, 
at this density the vacuum as connected medium also still exists (see Fig. E^). 




To prevent infinite connecting vacuum clusters, a much higher hadron density 
is needed. Measured in hadronic size units, the vacuum disappears for 

1 24 

n c = — ~ 0.56 fnT 3 , (12) 
Vh 

schematically illustrated in Fig.[7|D. If we assume that at this point, the medium is 
of an ideal gas of all known hadrons and hadronic resonances, then we can calculate 
the temperature of the gas at the density n c : n rcs (T c ) = n c implies T c ~ 170 MeV, 
which agrees with the value of the deconfinement temperature found in lattice QCD 
for n = 0. 

We can thus use percolation to define the states of hadronic matter. At low 
density, we have a hadron gas, which at the percolation point n c turns into connected 
hadronic matter. When this becomes so dense that only isolated vacuum bubbles 
survive, at n C) it turns into a quark-gluon plasma. This approach provides the correct 
values both for the density of standard nuclear matter and for the deconfinement 
transition temperature. 

Such considerations may in fact well be of a more general nature than the prob- 
lem of states and transitions in strong interaction physics. The question of whether 
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symmetry or connectivity (cluster formation) determines the different states of 
many-body systems has intrigued theorists in statistical physics for quite a long 
time^-4 The lesson learned from spin systems appears to be that cluster formation 
and the associated critical behaviour are the more general feature, which under 
certain conditions can also lead to thermal criticality, i.e., singular behaviour of the 
partition function. 



References 

1. I. Ya. Pomeranchuk, Doklady Akad. Nauk. SSSR 78 

2. R. Hagedorn, Nuovo Cim. Suppl. 3 (1965) 147; 
Nuovo Cim. 56A (1968) 1027. 

3. K. Bardakci and S. Mandelstam, Phys. Rev. 184 (1969) 1640; 
S. Fubini and G. Veneziano, Nuovo Cim. 64A (1969) 811. 

4. N. Cabbibbo and G. Parisi, Phys. Lett. 59 B (1975) 67. 

5. H. Satz, Fortsch. Physik 33 (1985) 259. 

6. M. Asakawa and T. Hatsuda, Nucl. Phys. A 610 (1996) 470c. 

7. K. Wilson, Phys. Rev. D10 (1974) 2445. 

8. M. Creutz, Phys. Rev. D 21 (1980) 2308. 

9. F. Karsch and E. Laermann, in Quark-Gluon Plasma III, R./ C./ Hwa (Ed.), World 
Scientific, Singapore 2004. 

10. L. D. McLerran and B. Svetitsky, Phys. Lett. 98 B (1981) 195 and Phys. Rev. D 24 
(1981) 450. 

11. J. Kuti, J. Polonyi and K. Szlachanyi, Phys. Lett. 98B (1981) 199. 

12. B. Svetitsky and L. G. Yaffe, Nucl. Phys. B 210 [FS6] (1982) 423. 

13. F. Karsch and E. Laermann, Phys. Rev. D 50 (1994) 6954. 

14. F. Karsch, E. Laermann and A. Peikert, Phys. Lett. B 478 (2000) 447; 
F. Karsch, Lect. Notes Phys. 583 (2002) 209. 

15. A. Peikert, F. Karsch and E. Laermann, Nucl. Phys. B 83-84 (2000) 390. 

16. M. Asakawa and T. Hatsuda, Nucl. Phys. A 610 (1996) 470c. 

17. J. Engels et al., Z. Phys. C 42 (1989) 341. 

18. E. Shuryak and I. Zahed, Phys. Rev. C 70 (2004) 021901. 

19. J. Engels et al., Z. Phys. C 42 (1989) 341. 

20. V. Goloviznin and H. Satz, Z. Phys. C 57 (1993) 671. 

21. F. Karsch, A. Patkos and P. Petreczky, Phys. Lett. B 401 (1997) 69. 

22. M. A. Halasz et al., Phys. Rev. D 58 (1998) 096007. 

23. Z. Fodor and S. Katz, JHEP 0203 (2002) 014. 

24. P. de Forcrand and O. Philipsen, Nucl. Phys. B 642 (2002) 290. 

25. M.-P. Lombardo, Phys. Rev. D 67 (2003) 014505. 

26. C. R. Allton et al., Phys. Rev. D 68 (2003) 014507. 

27. M. E. Fisher, Physica 3 (1967) 255. 

28. H. Satz, Nucl. Phys. A 642 (1998) 130c. 

29. V. Magas and H. Satz, Europ. Phys. J. C 32 (2003) 115. 



